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Abstract: We present an improved determination of the CP and CPT violation param- 
eters Re(e) and Im(5) based on the unitarity condition (Bell-Steinberger relation) and on 
recent results from the KLOE experiment. We find Re(e) = (159.6 ± 1.3) x 10~ 5 and 
Im(<5) = (0.4 ± 2.1) x 10 -5 , consistent with no CPT violation. 
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1. Introduction 

The three discrete symmetries of quantum mechanics, charge conjugation (C), parity (P) 
and time reversal (T), are known to be violated in nature, both singly and in bilinear com- 
binations. Only CPT appears to be an exact symmetry of nature. Exact CPT invariance 
holds in quantum field theory, which assumes Lorentz invariance (flat space), locality and 
unitarity [1]. Testing the validity of CPT invariance therefore probes the most fundamen- 
tal assumptions of our present understanding of particles and their interactions. These 
hypotheses are likely to be violated at very high energy scales, where quantum effects of 
the gravitational interaction cannot be ignored [2]. On the other hand, since we still lack a 
consistent theory of quantum gravity, it is hard to predict at which level violation of CPT 
invariance might become experimentally observable. 

The neutral kaon system offers unique possibilities for the study of CPT invariance. 
From the requirement of unitarity, Bell and Steinberger have derived a relation, the so- 
called Bell-Steinberger relation (BSR) [3]. The BSR relates a possible violation of CPT 
invariance (m K0 7^ m-^ and/or r^ 7^ ^~ko) in the time-evolution of the K°-K° system 
to the observable CP-violating interference of Kl and K$ decays into the same final state 
/. Strictly speaking, evidence of CPT violation found via the BSR could just be a failure 
of the unitarity assumption. However, unitarity is also one of the main hypotheses of 
the CPT theorem; thus the BSR allows a test of the basic assumptions of quantum field 
theories. 
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In this work we use recent results from the KLOE experiment to improve the deter- 
mination of the phenomenological CP- and CPT-violating parameters Re(e) and Im(<5) by 
means of the BSR. Our analysis benefits in particular from three new measurements: i) 
the branching ratio BR^i^— >7r + 7r~) [4], which is relevant to the determination of Re(e); ii) 
the new upper limit on BK(Ks— >7r°7r 7r°) [5], which is necessary to improve the accuracy 
on Im(<5); and iii) the measurement of the semileptonic charge asymmetry As [6], which 
allows, for the first time, the complete determination of the contribution from semileptonic 
decay channels without assuming unitarity. 

A determination of Re(e) and Im(5) using the BSR was performed by CPLEAR [7] 
in 1999. In the analysis of Ref. 7, some of the parameters of the semileptonic channels 
were evaluated together with Re(e) and Im(<5) from a combined fit to the time-dependent 
semileptonic asymmetries imposing the constraint of the BSR. A recent update of the 
determination of Re(e) and Im(<5) is given in Ref. 8. In this latter analysis, however, some 
of the results of the CPLEAR fit [7], which used the unitarity constraint, have been used 
as input again to the BSR. It is not clear whether this fact was accounted for in [8]. 

Our presentation is organized as follows. In section |2| we outline the meaning of the 
BSR and the approach used to maximize the sensitivity obtainable from the presently 
available data. In section |3] we examine experimental inputs, including a re-analysis of 
the semileptonic channels, and obtain a best estimate for errors and correlations. The 
extraction of Re(e) and Im(5) is discussed in section ^. 

2. Theoretical framework 

Within the Wigner-Weisskopf approximation, the time evolution of the neutral kaon system 
is described by [9] 

i^¥(t) = H*(t) = (m- l -Tj vf(i) , (2.1) 

where M and E are 2x2 time- independent Hermitian matrices and ^f(t) is a two-component 
state vector in the K°-K° space. Denoting by my and E^ the elements of M and E in the 
basis, CPT invariance implies 

mn = m-22 (or m K0 = m^ ) and r n = E 2 2 (or T R0 = T^ ) . (2.2) 

The eigenstates of eq. ( |2.1[ ) can be written as 

Ks ' L = /o n m ft 1 + es ^ K ° ± ( 1 ~ *S,l) K°) , (2.3) 

V 2 ! 1 + \ e S,L\ 2 ) 

-ilm(m 12 ) ~ flm(ri 2 ) ± \ (m^ - m K o - \ (T^ - E^o)) 
6S ' L m L -m s + i(T s - T L )/2 

= e±6, (2.4) 

such that 5 = in the limit of exact CPT invariance. 

Unitarity allows us to express the four elements of E in terms of appropriate combina- 
tions of the kaon decay amplitudes Af. 

r« =Y,Mf)Mf)*> h3 = 1,2 = K°,K°, (2.5) 
/ 
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where the sum runs over all the accessible final states. Using this decomposition in eq. (|2.4| ) 
leads to the BSR: a link between Re(e), lm(5), and the physical kaon decay amplitudes. In 
particular, without any expansion in the CPT-conserving parameters and neglecting only 
0(e) corrections to the coefficient of the CPT-violating parameter 8, we find 

(^i^+itan^sw) (J^-iIm(6^j = ^-1^ ^ ^(/)^*(/), (2.6) 

where (psw = arctan (2(m_L — Tns)/(Ts — IV)). We stress that, in contrast to similar 
expressions which can be found in the literature, eq. ( [2.6| ) is exact and phase-convention 
independent in the exact CPT limit: any evidence for a non-vanishing Im(<5) resulting from 
this relation can only be attributed to violations of: i) CPT invariance; ii) unitarity; iii) 
the time independence of M and T in eq. (|2.1| ). 

The advantage of the neutral kaon system is that only a few decay modes give signifi- 
cant contributions to the r.h.s. in eq. ([2.6|). Only the 7T7r(7), tttttt and irlv modes turn out 
to be relevant up to the 10~ 7 level 1 . 

The products of the corresponding decay amplitudes are conveniently expressed in 
terms of the a, parameters defined below. 

2.1 Two-pion modes 

For two-pion states, we define the ratios «j as: 

a, = ±- (A L (i)A* s (i)) = m BR(K S - i) , i = ttV , 7r+7r-( 7 ), (2.7) 
I s 

where tt + tt~(^) denotes the inclusive sum over bremsstrahlung photons, and (...) indi- 
cates the appropriate phase-space integrals. The rji parameters in eq. fl2.7| ) are the usual 
amplitude ratios: rji = Ai{i) / 'As{i). 

The contributions from tt + tt~j direct-emission (DE) amplitudes not included in a n + n -^ 
are collected together in the term 

a TTTDE = a i"T7El-S + a T7r7 E1 „ L + Q 7r7T7 DExDE ; (2-8) 

where 

= ((^L(vr7r7)^5(7T7r7Ei)) + (^(^^1)^(^7))) 

= AB(K S -> 7T7T7DE) V+- + (J1+-J ~ V+-) BR(K S -> VT7T7). (2.9) 

Here Al si^^l) an d ^L,s (^Tei) denote the leading bremsstrahlung and the electric- 
dipole DE amplitudes, respectively. Their interference cannot be trivially neglected. RR{K$ 
— > ttttj) indicates the branching fraction for decays with the emission of a real photon 



Note that all quantum numbers of a chosen final state / have to be equal between Ks and Kl decays 



in order to allow for interference between the two amplitudes in the r.h.s. of eq. (2.6) 
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with minimum energy equal to the cut used in the corresponding rj-\ ~ measurement. 

AB(K S -» 7rvr7 DE ) = BR(K S -» vr7T7) cxp - BR(K S -* irir I ) & ~ m is the DE contribution to 
the BR, obtained subtracting the computed bremsstrahlung spectrum from the measured 
spectrum (see appendix |A]). 

We have generically denoted by 07T7T7 DExDE the contribution arising from the product 
of two DE amplitudes (electric or magnetic). Given the strong experimental suppression 
of DE amplitudes, this term turns out to be safely negligible being of 0(1O -8 ) or less [11]. 

2.2 Three-pion modes 

For three-pion states we define 

Oi = -^(A L (i)A s (i)} = —ri* BR(K L - i), i = 3vr° , A+tt"^. (2.10) 
L S t Kl 

Note that in this case the amplitudes are not necessarily constant over the phase space. 
As a result, the rji appearing in eq. ( [2.10 ) should be interpreted as appropriate Dalitz-plot 
averages. In particular, the ir + Tr~ir° final state is not a CP eigenstate. For decays to 
7T + 7r~7r° the rj parameter can be expressed as 

_ (A* l {tt + tt-tt°, CP = +l)^ 5 (vr + 7r-vr , CP = +1)> 
V+ ~° ~ (|^ L (7r+vr-vrO)| 2 ) + 

(AliTT+TT-TT , CP = -l)As(7T + 7T-7r , CP = -1)) 
(\A L (7T+7T~^)\ 2 ) 

The experimental bounds on rj^ o reported by CPLEAR [10] correspond to this average, 

when neglecting the contribution of Al{^ tt° ', CP = +1). This is indeed a good ap- 
proximation, because the Kl decay amplitude to a CP = +1 7r + 7r~7r° state is suppressed 
both by CP violation and a centrifugal barrier. Given the poor direct experimental infor- 
mation on 77000 ) i n the neutral case it turns out to be more convenient to set a bound on 
1 0^.0^0^.0 1 by means of the relation 

ItVW)! 2 = — BR(K L -> 3^°) x BR(K S -» 3vr°) . (2.12) 

This relation is based on the well-justified assumption that the Kl (Kg) decays to 3ir° are 
dominated by a single CP conserving (violating) amplitude with the same behaviour over 
phase space [11]. 

2.3 Semileptonic modes 

In the case of semileptonic channels, the standard decomposition is [12] 

A{K° -► %+rv) = A* (l + y*)(x+ - x_)* , 
A(K°^n + ru) = A* {l + y*) , 

A(K° -» n'l + u) = Ao(l - y)(x+ + z_) , (2.13) 
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where x+ (%-) describes the violation of the AS = AQ rule in CPT conserving (violating) 
decay amplitudes, and y parametrizes CPT violation for AS = AQ transitions. Assuming 
lepton universality, and expanding to the first non-trivial order in the small CP- and 
CPT-violating parameters, one obtains 

Y^iMl{kIv)A* s {i:Iv)) = 2T{K L -> ir£v) (Re(e) - Re(y) - i(lm(x + ) + Im(£))) 

= 2T(K L -> irtv) ((A s + A L )/4 - i(Im(x+) + Im(<5))) . (2.14) 

The dependence of Re(y) has been eliminated by taking advantage of the relation Re(e) — 
Re(y) = (As + A£) / '4 [12], where Al,s are the observable semileptonic charge asymmetries. 
The parameter Im(x + ) can be measured using the appropriate time-dependent decay dis- 
tributions [13], while Im(<5) is one of the two outputs of the BSR. In order to get rid of the 
explicit Im(5) dependence, it is convenient to define 



— V (A L (nlv)A s (irlv)) + 2i^ BR(K L -» 7rft/)Im(<J) 



s 7Tu TR i 

= 2^BR(K L ^ir£u)((A s + A L )/4-ilm(x + )) . (2.15) 
T K L 

2.4 Determination of Re(e) and Im(<5) 

The on parameters defined in eqs. (|2.7| ), ( ^q) , ( |2.10|) , and ( 2.15 ) can be determined (or 
bounded) in terms of measurable quantities. Taking into account these definitions (in 
particular the non-standard expression for a n £ u ), the solution to the unitarity relation in 
eq. ( |2.6P is 

Re(e) 

1 17 ^ " ^ )(^ e ^M (2 . 16) 

;i-K)tan0sw + / \ E i Im (ai)/ 




= BR(K L -» vr£i/), and 
iV = (1 + k) 2 + (1 - k) 2 tan 2 </>sw - 2 6 k(1 + «) . (2.17) 



3. Experimental input to the a parameters 

The experimental inputs needed for the determination of the oti are the Kl and Ks branch- 
ing ratios, the amplitude ratios rji, and the Kl and Ks lifetimes. All experimental inputs 
used in the determination of the decay amplitudes are summarized in table |]. For the deter- 
mination of the neutral kaon decay rates and the Kl lifetime we combine the measurements 
listed below. 

1. The absolute K L BR's from KLOE [21]; 

2. The KLOE K L lifetime [20]; 

3. The KLOE ratio BR(K L -> tt + tt-)/BR(K l -> irfiu) [4]. This result is inclusive 
of final state radiation, therefore the DE contribution in the process Kl^tt + tt^ is 
subtracted using the result from Ref. [18]; 
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Table 1: Input values to the Bell-Steinberger relation. Results from [14] are those evaluated 
without assuming CPT invariance. The KLOE average and the average are described in Sec. 3 
and Sec. 3.3, respectively. 



4. The precise KLOE determination of R s = BR(K S -» 7r+7r _ )/BR(As -» vr 7r°) [16], 
and the value of Rs /Rl from the world average of Re(e' /e) [14] , with Rl = BR(i^L — > 
7t + tc~)/BR(Kl —i- 7T°7T°). These measurements are used to determine the value of 
BR(i^ L ^7r 7r°). 



The combination of all of the above listed measurements, referred to as the KLOE aver- 
age, accounts for all correlation effects, and is obtained by renormalizing the sum of Kl 
branching ratios to 1 — BR(i£^ — > 77). This procedure yields very small corrections to the 
published BR values [21]. The results with errors and correlation coefficients are given in 
table [| 

For the K$ lifetime we use the average tk s = 0.08958 ± 0.00006 ns [14] obtained 
without assuming CPT invariance. 
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Value 


Correlation coefficients 


BR(K L - 




0.4009(15) 


1 






BR(K L - 




0.2700(14) 
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1 




BR{K L - 


>7r + 7T _ 7T ) 


0.1262(11) 
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1 
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-0.15 
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-0.06 -0.23 1 


BR(K L - 
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8.48(10) xl0~ 4 


-0.14 


0.49 


-0.06 -0.23 0.97 1 


tk l (ns) 




50.84(23) 


0.16 


0.21 


-0.26 -0.13 0.07 0.06 1 


Rs 




2.2549(54) 


0.00 


0.00 


0.00 0.00 0.00 -0.21 0.00 1 



Table 2: Values, errors, and correlation coefficients for all parameters included in the KLOE 
average: the dominant Kl BR's, Kl lifetime and Rs- 



3.1 Two-pion modes 

The 7T7r(7) terms are the largest. The value of a nn of eq . (|2 . 7f) is evaluated as: 



1/2 

a mt =[-^-x BR(K L -> tttt) x BR(K s -> tttt 



TR S v F!R| IV. , RRi TV.. \ •>■■'■■■> 

The values of BR(Ks—>tt + tt~) and BR(K$— >7r 7r°), given in table 0, are determined from 
R s using the constraint BR{K S -> vr+vr") + BR(K S -> 7r°vr ) + BR(i^ s -> vr£i/) = 1; the 
value of BR(i^5 — ► 7r£^) is determined from the KLOE measurement of BR(i£"5 — * Trev) [6] 
assuming lepton universality. For BR(Kl — > 7r + 7r~) and BR(i^i — > 7r°7r°), we use the 
values from the KLOE average, which are in agreement with recent measurements from 

KTeV [15]. Finally, the values of ^ and 0oo> the phases of r/^ and r/oo, are taken from 

the PDG fit [14] without assuming CPT invariance. 

Figure [j] shows the 68% and the 95% CL contours for a n + w - and a n o n o. We find: 

a7r+7r - = ((1.115 ± 0.015) +*(1.055± 0.015)) x 10" 3 , 
a n07T o = ((0.489 ± 0.007) + i(0.468 ± 0.007)) x 10~ 3 . 



Following the discussion in section 2T, we determine a 7r + 7r - 7DE from the measurement 

of 77-1 7 from Ref. [17] and BR(Ks — > 7r + 7r~7) from the measurement of Rs,-y = BR(Ks — > 

7r + 7r~ 7) /BR(K s — > vr + 7r _ ) [18]. The photon energy threshold was 20 MeV for both 
measurements. The value of R^~ m = (BR(K S -> ttttj)/BR(K s -> 7rvr)) th " IB = (0.700 ± 
0.001) x 10~ 2 is calculated using the KLOE MC generator, which is described in Ref. 19. 
Figure ^(a) shows the 68% and the 95% CL contours for a w + 7r -^ DE . We find a very small 
contribution: 

a7T+n - 7DE = ((5 ± 7) + i(6 ± 7)) x 10- 7 . 
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Figure 1: Representation of a^^- and a^a^o in the complex plane. In each case, the two ellipses 
represent the 68% and the 95% CL contours. 
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Figure 2: Bounds of the a values for three-body decays of K$, K^. Note that the same scale is 
used for all plots. 

3.2 Three-pion modes 

Until the recent results from KLOE [5] and NA48 [8] became available, the limiting con- 
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tribution to the CPT test was due to the 3ir° final state. The upper limit on lo^o^o^o | is 
evaluated from eq. ( fEl2l) using the KLOE upper limit BR(i^5^3vr°)< 1.5 x 10~ 7 at 95% 
CL [5] and the value of BR(1^^37r ) from the KLOE average. The phase <^ooo is taken as 

uniform in {0, 2ir}. For a n + 7T - 7T o, eq. ( |2.10| ), we use Re(r/_| o) and Im^ o) a s measured 

by CPLEAR [10] and the value of BR(i^i^7r + 7r _ 7r ) from the KLOE average, all given in 
table 0. 

The 68% and the 95% CL contours for a n + n -^o, and the 95% CL contour for a^o^o^o 
are shown in figures ||(b) and §(c), respectively. We find: 



o 



((0 ± 2) + i(0 ± 2)) x 10~ 6 , 
6 at 95% CL . 



|a 7r o 7r o 7r o I < 7 x 10 
3.3 Semileptonic modes 

For the determination of a n e u from eq. (2.15), we combine the KLOE measurement [6] of 
the Ks semileptonic charge asymmetry As for K$ semileptonic decays, the PDG aver- 
age [14] for Kl semileptonic charge asymmetry Al, and the CPLEAR time-dependent 
measurement of iT and K° semileptonic rates [13]. The original CPLEAR result is given 
in table |3[ 





value 


Correlation coefficients 


Re{6) 


(3.0 ± 3.4) x 10" 4 


1 


lm(S) 


(-1.5 ±2.3) x 10~ 2 


0.44 1 


Re(x_) 


(0.2 ± 1.3) x 10~ 2 


-0.56 -0.97 1 


Im(x + ) 


(1.2 ±2.2) x 10~ 2 


-0.60 -0.91 0.96 1 



Table 3: Values, errors, and correlation coefficients for Re(<5), Im(<5), Re(x_), and Im(ir + ) measured 
by CPLEAR. 



We have improved this result by adding the measurement of As — Al = 4[Re(<5) + 
Re(x_)] = (-2 ± 10) x 10" 3 [6,14]. The results, referred to as the K a average, are given 
in table |I[ Finally, for BK(Kl — > nlv) we use the sum of K e ^ and branching ratios 





value 


Correlation 


coefficients 


Re(5) 


(3.4 ±2.8) x 10~ 4 


1 




lm(S) 


(-1.0 ±0.7) x 10~ 2 


-0.27 1 




Re(x_) 


(-0.07 ±0.25) x 10~ 2 


-0.23 -0.58 


1 


Im(x + ) 


(0.8 ± 0.7) x 10~ 2 


-0.35 -0.12 


0.57 1 


A S + A L 


(0.5 ± 1.0) x 10~ 2 


-0.12 -0.62 


0.99 0.54 1 



Table 4: Values, errors, and correlation coefficients for Re(<5), Im(i5), Re(x_), Im(a;+), and A$ + Al 
obtained from a combined fit (K 13 average). 



from the KLOE average. Figure [|(d) shows the 68% and the 95% CL contours for a n £ v . 
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We find: 



a n£u = ((0.3 ± 0.6) + *(— 1-8 ± 1.8)) x 10~ 5 



4. Results 

Inserting the values of the a parameters into eq. (2.16), we obtain 2 : 



Re(e) = (159.6 ± 1.3) x 10 



-5 



Im(6) = (0.4 ±2.1) x 10" 



(4.1) 



where all correlations among the input data are taken into account, including that from 
the direct determination of Im(<5) from the semileptonic decays. The complete information 
is given in table ||. 





value 


Correlation coefficients 


Re(e) 


(159.6 ± 1.3) x 10" 5 


1 


lm(S) 


(0.4 ±2.1) x 10~ 5 


-0.17 1 


Re(5) 


(2.3 ± 2.7) x 10" 4 


0.20 -0.22 1 


Re(x_) 


(-2.9 ±2.0) x 10~ 3 


-0.25 0.37 -0.49 1 



Table 5: Summary of results: values, errors, and correlation coefficients for Re(e), Im(<5), Re(<5), 
and Re(x_). 



The allowed region in the Re(e), Im(<5) plane at 68% CL and 95% CL is shown in the 
left panel of figure H. A small correlation between Re(e) and Im(<5) is evident, which is 



0.5 



-0.5 



Imd (10 4 ) 



95% CL 
68% CL 




Re£ (10 2 ) 



0.155 



0.160 



0.165 



10 



•10 



AT (10 18 GeV) 



95% CL 
68% CL 



18 



AM (10 GeV) 



-10 



10 



Figure 3: Left: allowed region at 68% and 95% CL in the Re(e), lm(6) plane. Right: allowed 
region at 68% and 95% CL in the AM, Ar plane. 

due to the semileptonic term. With the new KLOE data used in the present analysis, the 
process giving the largest contribution to the size of the allowed region is now Kl — > 7r + 7r~, 

2 The accuracy on Re(e) improves by about 30% with inclusion of the KTeV measurement of BR(A~l — > 

TT+O [15]. 
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through the uncertainty on cp^ Besides the tttt final states, only the semileptonic term 

gives an appreciable contribution (~ 10%) to the error on Im(<5). Our results, Eq. 4.1, are 
a significant improvement over those of CPLEAR [7]: 



Re(e) = (164.9 ± 2.5) x 10~ 5 , Im(S) = (2.4 ± 5.0) x 10~ 5 . 

Note also that the central value of Re(e) is quite different. This is due to the new mea- 
surement of BR(.K£ — >7T + 7T~) [4]. 

The limits on Im(5) and Re(5) can be used to constrain the mass and width difference 
between K° and K° via 

6 = i{m K0 -m^) + \{T K0 -T^) ^ {1 + Q{e)] 

The allowed region in the AM = (m K — m^ ) , Ar = (r^ — ) plane is shown in the 
right panel of figure |3[ The strong correlation reflects the high precision of lm(5) compared 
to Re(<5). Since the total decay widths are dominated by long-distance dynamics, in models 
where CPT invariance is a pure short-distance phenomenon, it is useful to consider the limit 
r^ = r^ - In this limit (i.e. neglecting CPT- violating effects in the decay amplitudes), 
we obtain the following bound on the neutral kaon mass difference: 

-5.3 x 10~ 19 GeV < m R0 - m F < 6.3 x 10~ 19 GeV at 95 % CL . 

Our result represents a significant improvement with respect to \m K0 — m-j^ \ < 12.7 x 
10~ 19 GeV at 90% CL, obtained by CPLEAR [10]. 
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A. The 7T7T7 contribution to the ctj 

By construction, the leading contribution of the ttttj state to the unitarity sum, namely 

the interference of the Kl and K$ bremsstrahlung amplitudes, is included in ( 7 ). The 

largest sub-leading term missing in ( 7 ) is the DE-bremsstrahlung interference, which 

we include in aw7 E1 _ s(L) of eq. (p.8[) . To evaluate this contribution, we introduce the total 
(IB+DE) amplitude ratio 



\ ^l(vT7T7ib+Ei) I / / \ /a -i \ 

^l s (7r7r7iB+Ei) 



This ratio is an observable quantity which can be measured in an interference experi- 
ment [17]. As explicitly indicated, r/^ 7 depends on the photon energy and it can be 
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decomposed into the energy-independent parameter ry_| and the direct-CP-violating term 

T In the Ej — ► limit e+_y ac [11]. Using eq. ( |A.1| ), we can write 



r 5 



— ( r? + _ (As(^nOAs( 7r7r 7Ei) + ^(^^1)^5(^7)) 



+ 



+e' + _ 1 As(i™i)A* s (mti){\ + 5 E1 )) 
. AB{K S -> 7rvr7 DE ) + ^ J dE 1 e\ 



+-7 



dT(Ks — > iriT'y) 
~~dEZ 



(A.2) 



where <5ei = ^4s(7T7r7Ei)/^4s( 7r7r 7) is a very small overall correction factor which can be 
safely neglected and AB(K S -> 7rvr7 DE ) = BR(K S -> 7rvr7) ex P - BR(JT S -> 7rvr7) th - IB 
is the deviation of the observed K$ —* ttttj decay rate from that inferred from a pure 
bremsstrahlung spectrum. By construction, the integral on the right-hand side of eq. ( |A,2| ) 

is infrared safe. Since at present there is no evidence for a non- vanishing e' ^ ~ [17], in 

eq. ( |2.9| ) we replace this integral with the product (77-1 7 — rj^ ) x BR(i^5 — > irirj), where 

BR(i^5 — > ttttj) indicates the branching fraction for a real photon emission with minimum 
photon-energy cut equivalent to that used in the corresponding r/^ ~ measurement (E^ ut = 



20 MeV in Ref. 17). The contribution to eq. (|A.2[) generated by -Kl,s — > ^Try amplitudes 

and can be safely neglected. 



with Ery < E™ vanishes in the limit E^ nt 
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